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Abstract 



We reply to a critique of our paper arXiv:1210.5501 by the DAMA collaboration which appeared 



in 



arXiv:1210.6199 Our original claim that the observed background levels are likely to require a 
large modulation fraction of any putative signal holds. In fact, in light of DAMA's recent comment 
our claim is further corroborated. We identify the source of the discrepancy between our own 
analysis and DAMA's claimed levels of unmodulated background. Our analysis indicates that the 
background in the signal region as reported by DAMA is likely underestimated. 
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In a recent publication pQ we pointed out that the empirical verification of a certain special 
decay branch of 40 K remains outstanding and discussed the importance of such a measure- 
ment. In addition, we discussed the general role 40 K plays as an important background in 
the DAMA experiment [2|. Shortly after our paper appeared, the DAMA collaboration has 
criticized some of its findings [3]. There are three points of substance in their critique: 

Critique 1: The collaboration states that the contribution of this special decay of 40 K to 
the total 40 K contribution is only about 10%. They therefore criticize our discussion 
of this branch as "captious". 

Critique 2: The collaboration states that the content of nat K has been measured by in- 
vestigation of the double coincidence, is independent of the branching ratio into the 
ground state, and its average value over all crystals is 13 ppb. Statements based on 
20 ppb are therefore wrong. 

Critique 3: Finally, the collaboration claims an upper limit of 0.25 cpd/kg/keV for the 
unmodulated part of the signal in the 2-4 keV energy bin. Given a residual rate of 
0.0194 ± 0.0022 cpd/kg/keV [2\, this would allow for a modulation fraction of 6-10% 
which can be accommodated with many models of Dark Matter j4], contrary to our 
conclusions. 

The first two issues are rather trivial. The third issue is very important and its elucida- 
tion requires understanding the details, which we give in the appendix. We will show that 
the upper limit of So < 0.25 cpd/kg/keV is likely the result of an underestimation of the 
background in the signal region. 

Reply 1: Independently of DAMA, the measurement of the 40 K decay directly to the ground 
state of 40 Ar is a missing piece in the experimental verification of leading-order predic- 
tions of weak decays, being the only such EC decay known of its kind (third unique- 
forbidden). In [1] we calculated the theoretical prediction for the branching ratio and 
motivated its measurement. The total decay of 40 K is a serious background in the 
DAMA experiment and it was in fact one of our findings that this special branch is 
subdominant to the one where 40 K decays into the 40 Ar excited state, followed by a 
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1.46 MeV gamma-ray which escapes detection. In their comment [3] DAMA quotes 
(for the first time in print) that the former contribution is 10% of the total 40 K con- 
tribution. This is indeed in good agreement with what we found using the theoretical 
branching ratio and using the independent Monte Carlo simulation of Ref. [5]. This 
agreement gives further credence to our treatment of potassium decays in DAMA. 

Reply 2: Our results are presented for a whole range of nat K contaminations between 
1-100 ppb including the DAMA quoted upper limit of 20 ppb, which appeared in 
the official publication describing the DAMA apparatus [B]. The collaboration now 
provides an average contamination level of 13 ppb in their comment j3]. We could not 
find this number in any of their published works but it can be found on slide number 8 
in a talk by Nozzoli [7J . More importantly, no data is presented to support this number, 
nor does the collaboration provide the uncertainty associated with it. Nevertheless us- 
ing 13 ppb we find that the required signal modulation fraction is above 20%. In [T] we 
identified a way to measure the contamination level in a manner that is free from the 
Monte Carlo modeling required by the coincidence method. Considering the immense 
importance of the precise amount of potassium background to any interpretation in 
terms of a dark matter signal the lack of details so far provided by the collaboration 
is unsettling. 

Reply 3: Again, we could not find any mention of an upper limit of 0.25 cpd/kg/keV for 
the unmodulated signal part in any of DAMA's published works. We therefore have to 
resort to slide 8 of a talk by Nozzoli [7J where the limit is quoted and a curve for the 
background model is reported. In their comment, DAMA says that they account for 
this limit in [I] when assessing the viability of various DM scenarios. In Fig. [I] we show 
our reproduction of this background, which is a Gaussian centered at 3.2 keV together 
with a linear fit as explained in the Appendix. As is clear from the plot, the linearly 
rising function is not supported by the data above 4 keV (the statistical error bars are 
shown but are barely visible.) More importantly, the decrease towards lower energies 
strongly supports the notion that the background present in the signal region 2-4 keV 
is underestimated. As argued in our original work [1J, a flat background component 
is expected because it is a universal feature of (3~ decays for small electron velocities. 
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Indeed, adding a flat component greatly improves the fit as shown in Fig. [Tj This results 
in a significantly larger background estimate in the signal region (see the Appendix for 
details). In light of this alternative estimate, the required modulation fractions quoted 
in our work [T], where a relatively low constant value of 0.85 cpd/kg/day was used, 
are in fact seen to be conservative. 

What should be clear from this reply is that a more detailed and careful discussion of 
backgrounds by the DAMA collaboration is called for. That was among the main conclusions 
in pQ where we also identified measures which could help to clarify some of the uncertainties. 
We invite the interested readers to consult the appendix where more details are presented 
and to reproduce our results themselves. 
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FIG. 1. The dots (with error-bars) are the single-hit rate reported by DAM A in [6]. The thick 
red curve shows our reproduction of the background curve presented by DAMA in Ref. [7j . The 
dashed curve from 6.2 keV to 10 keV is the continuation of the linear trend. The linear trend is 
obtained by fitting the data between 5.3 keV and 10 keV and results in excellent agreement with 
the reported background [7j, as explained in the Appendix. Using the same Gaussian, but with 
a better model for the background above 5 keV (flat + linear), we obtain the solid green curve. 
In blue, we show the resulting curve in the case when nat K contamination level is increased to 14 
ppb instead. We emphasize that we did not fit to the data below 5 keV, but used the Gaussian as 
reported by the collaboration in [7] . 
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Appendix A: Details of the analysis 



The unmodulated background presented on slide 8 of Nozzoli [7] is reproduced accurately 
by fitting a linear function to the data points between 5.3 keV and 10 keV and adding to it 
a single Gaussian whose parameters are taken directly from ref. [7j, 



A ( (E-E 



\2 



BKG Nozzo ii(£) = ^ exp ( - % ) + slope x E + intercept. (A-l) 
v27r<7 2 \ ^ a J 

The amplitude of the Gaussian is A = 0.64 which is in good agreement with the results 
of [5] with a nat K contamination of 13 ppb. The center and spread of the Gaussian are 
in good agreement with the expected 3.2 keV energy deposition (E KiQ = 3.15 keV) and 
resolution (er = 0.618 keV). One would naively expect that the quoted upper limit So < 
0.25 cpd/kg/keV was derived by subtracting this simple background model from the data. 
However, in contradiction to what the DAMA collaboration is quoting this yields only Sq ~ 
0.14 cpd/kg/keV in the 2-4 keV energy bin as can be easily checked. 



Given the negligible statistical error bars, it is seen in Fig. [I] that the linear fit in Eq. (A-l) 



provides an inadequate description of the data between 5-10 keV where little DM signal is 
expected. A flat background up to about 7 keV where it is broken and followed by a linearly 
rising one provides a much better agreement with the dataQ 

0.866 E < 7.05 keV 
BKG linear+flat (£) = ^ (A-2) 

0.491 + 0.053E E > 7.05 keV 

Such a model is also physically well-motivated. The background from /^"-emitters is entirely 
flat at low energies. The rise above 7 keV is more difficult to interpret because the DAMA 
collaboration has not released the the immediate spectrum above 10 keV. For example, it 
could be a broad Iodine escape peak between 14-18 keV originating from external 210 Pb con- 
tamination. Though such features have been measured previously with Nal(Tl) crystals [S], 
without further insight into the spectrum this remains a speculation. 

1 We have also tried a more general model involving two linear fits, but the results are qualitatively 
unchanged - the fit between 5-7 keV prefers to be flat. 
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FIG. 2. The solid-red curve shows the background model presented in ref. [7]. The dashed-red curve 



represents a full fit to the unmodulated DAMA data composed of a Gaussian component, a flat 
component, and a linearly rising component. The center of the Gaussian, its spread, the flat piece, 
and the slope and intercept of the linear function were allowed to float (total of 5 parameters). We 
stress that the amplitude of the Gauassian was fixed to a level which corresponds to a 13 ppb of 
nat K contamination as quoted by the collaboration. 

We have also verified that a more general fit to the entire data (including the width and 
center of the Gaussian) does not result in any substantial change to the flat component 
which becomes 0.86 cpd/kg/keV, in agreement with what we used in our original work pp. 
In fact, this model results in a disturbingly good fit to the entire data set as shown in Fig. [2] 
We emphasize that this fit is done with a fixed amplitude for the Gaussian determined by a 
contamination level of nat K of 13 ppb as quoted by the collaboration. Only five parameters 
were allowed to float in the fit: the centre of the Gaussian; the resolution; the flat component; 
and the slope and intercept of linear trend. This surprising result shows how volatile the DM 
interpretation is to small changes in the background model and it again emphasizes the need 
for a more thorough discussion of backgrounds by the DAMA collaboration. In generating 
Fig. [I] above we attempted to remain as accommodating as possible and used the Gaussian 



is the extrapolation of the linear function as explained in the caption of Fig. [TJ The solid-blue line 
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curve quoted in [7j for modeling the bump and not the fit shown in Fig. |2} 

We notice that the residual signal that would result from Fig. [T] is so small that it would 
require a modulation fraction in excess of 50%. This is higher than the more conservative 
value we quoted in our paper P of a little over 20%. The reason for this is that the Gaussian 
model we assumed in our paper was based on the resolution as given by DAMA [6], which 
tends to be a little broader than what is actually observed in the bump. This results in 
slightly less potassium background in the signal region of 2-4 keV. Therefore a larger residual 
is possible and the modulation fraction is correspondingly smaller. 

Finally, we stress that the purpose of the background analysis we present is not to claim a 
full understanding of the backgrounds. Our goal, as it was in our paper, is to call attention 
to this important background and encourage the collaboration to present a full account of its 
detailed understanding. That being said, the remarkable agreement of our simple background 
model with the data must give pause to anyone who wishes to interpret the DAMA results 
in terms of a signal of dark matter. 
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